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Abstract 

We study theoretically the parametric down conversion and squeezing of microwaves using cavity 
quantum electrodynamics of a superconducting Cooper pair box (CPB) qubit located inside a 
transmission line resonator. The non-linear susceptibility \2 describing three-wave mixing can be 
tuned by dc gate voltage applied to the CPB and vanishes by symmetry at the charge degeneracy 
point. We show that the coherent coupling of different cavity modes through the qubit can generate 
a squeezed state. Based on parameters realized in recent successful circuit QED experiments, 
squeezing of 95% ~ 13dB below the vacuum noise level should be readily achievable. 
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Squeezed states are a valuable tool to usefully manipulate the Heisenberg uncertainty 
principle by reducing the quantum fluctuations of a certain variable of interest at the expense 
of increased uncertainty for its conjugate variable. Using squeezed states one can perform 
very quiet measurements much below the vacuum noise level [lj. Squeezed states manifest 
the quantum coherent nature of light and provide a chance to beat the standard quantum 
limit by preferentially doing an experiment using the squeezed quadrature alone 
Squeezed states have been experimentally observed in a nonlinear optical cavity experiment 
Q Q. Re C en tly , t he tey of „ g in a M gh - Q cav lty « considered Q. Upon t ne 
injection of a high-energy photon, a nonlinear optical medium can coherently generate two 
photons, the sum of whose frequencies is equal to that of the high-energy photon via optical 
parametric down conversion (PDC). If one injects low-energy photons instead, one may 
induce second harmonic generation, which also forms a squeezed state. In addition to this 
three-wave mixing, four-wave mixing can be used to generate squeezed states. In pioneering 
condensed matter experiments, the Josephson junction parametric amplifier was used in the 
microwave regime to produce (47 ± 8)% ~ 3dB squeezing below the vacuum level [3, 0, ^ 
via degenerate four wave mixing. There has been tremendous recent progress in realizing 
quantum optics physics in electrical circuits. It is now possible to experimentally reach the 
extreme strong coupling limit of cavity QED 11, [3| and to see very strong microwave 



non-linearities in high inductance small scale Josephson junctions whose Hamiltonian can 
be controlled with remarkable accuracy Q|. Coherent dynamics of a flux qubit coupled to 
a harmonic oscillator in a SQUID circuit has also been demonstrated l^j ]. 

Motivated by these experimental advances in 'circuit QED', we here study squeezing in 
a system consisting of a Cooper pair box (CPB) located inside a high Q coplanar waveg- 
uide resonator. Two different discrete photon modes (fundamental and first harmonic) are 
coupled through the CPB as shown in Fig. ^ and squeezing occurs via three- wave mixing. 
Compared to the cavity QED of atomic physics, circuit QED has the advantages of infinite 
transit time of the 'atom' (qubit) inside the cavity and much stronger coupling between 
qubit and photon. We emphasize that strong coupling means we only need a single 'atom'. 
The circuit QED system recently successfully demonstrated strong coupling (vacuum Rabi 
splitting) between a single photon and a qubit in all solid state system [ll|. While this archi- 
tecture offers a very quiet environment leading to excellent coherence of the qubit (7\ ~ 7/zs, 
T 2 * ~ 500ns) [l^, the solid-state environment still leads to qubit decay rates and dephasing 
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rates much larger than those in corresponding atomic physics microwave cavity QED exper- 
iments Hence it is crucial to investigate these environmental effects on the efficiency 
of squeezing. Using both numerical and analytical calculations based on the currently avail- 
able experimental parameters, we have estimated that the circuit QED system can readily 
produce about 95% squeezing, that is, 13dB below the vacuum noise level. 

We start with the following Hamiltonian to describe the coupled system of qubit and 
cavity photons for microwave circuit QED 
given by 



10, [111: H = H + E T . The Hamiltonian H is 



H = ~Y a z + hoJia\ai + huj 2 a 2 a 2 , (1) 

where E i = ^ Ej + E^, Ej is the Josephson coupling energy, E e \ = 4Ec*(l — 2N g ), E c 
the charging energy, N g the gate charge, and uj\ the angular frequency of the fundamental 
resonator mode and uj 2 = 2u\ is the first harmonic frequency. The coupling Hamiltonian Hi 
can be written 

Hi = - gi{a\ + a{) + g 2 {a\ + o 2 ) [1 - 2N g + sin 9 a x - cos 9 a z ) , (2) 

where gi represents the coupling strength between the qubit and the ith cavity photon mode, 
and 9 = tan~ l (Ej/E e \). At the charge degeneracy point (CDP), that is, N g = 1/2, Hi 
reduces to the Jaynes-Cummings Hamitonian, since cos 9 = —(4Ec/Ej)(l — 2N g ) = and 
sin# = 1. One can see that away from the CDP, we have couplings other than the Jaynes- 
Cummings term, whose strength linearly increases with the deviation. Because the only 
relevant cavity frequencies will be the fundamental and first harmonic, we expect degenerate 
PDC to occur in our system, where a single high-energy photon coherently generates two 
photons, each with half the frequency. We emphasize that this is degenerate PDC because 
the cavity only has discrete modes. In order to achieve squeezing via PDC in circuit QED, 
the system should be able to convert a single lv 2 photon into two uj\ photons, which requires 
a term of the form a|a|a2 in the effective Hamiltonian. This will result from the third order 
processes in terms of Hi, which will be in general negligible since gfg 2 /uf C 1 in typical 
experiments 11]. However we can resonantly enhance the process by tuning either u\ or uo 2 
close to Eqi- We have chosen the case uj 2 = E i- 

We first apply the following unitary transformation U\ to the Hamiltonian H 



U\ = exp 



X; -(a t -af)(l-2N g -cos9a z 



i=1.2 U ' 



(3) 
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This corresponds to shifting the centers of the harmonic oscillator coordinates X% = ai + a\ 
by {2gi/uji){l — 2N g — cos9<j z ). Subsequently, we apply the unitary transformation 



U 2 = exp 



gi sin 9 
2{E 01 - W) 



(at 



a — d\<7 



(4) 



Upon application of these two unitary transformations, we obtain after perturbative expan- 
sion in gi and g 2 , the following Hamiltonian: H e s = U 2 UiHUiU 2 = H + H' for E i = uj 2 , 

tti 1 ni t — -4-\ 9i sin 29 . 4. + _ 1, , s 

H = — -g 2 sin t){a 2 a + a 2 <7 ) H (a[a[a + a\a\a ). (b) 

2 2uj\ 

We define the energy detuning between the cavity photon frequency uj 2 and the qubit energy 
splitting Eq\ to be A = Eq\ —uj 2 and consider the case g 2 <C A -C uj\. It is this Hamiltonian 
which we will study numerically. However to develop an analytical understanding, we can 
apply the following additional unitary transformation f/3 

g 2 sin 9 



U 3 = exp 



{a^a — a 2 a + ^j 



(6) 



2A 

Finally, we obtain the following low-energy effective Hamiltonian through H e s = UsH e f[Ul, 

H cS = H + — [1 + a z (2n 2 + 1)] + -(a\a\a 2 + aia x a\)<j z , (7) 

where the second term on the right represents the Lamb and light shifts of the qubit splitting 
frequency and ( = (2gfg 2 sm9 sin29/Aui). The third term is the desired squeezing term. 
H c fi is exactly the standard Hamiltonian for degenerate optical PDC including the Lamb 
shift (3(. Note that the squeezing operator couples to the qubit state a z and the phase of the 
squeezed quadrature will shift by tt/2 if the qubit is placed in the excited state. 

In deriving the above result, we have neglected the effect of cavity damping 
J2i=i,2 \A* J^oo a «^L where b u represents the continuum modes outside the cavity and 
Ki the cavity loss rates. Unitary transformation of the damping terms leads to radiative 
atom damping and two-photon decay terms such as (y^C/^A) du aiaifo^. The 
two-photon decay rate is much smaller than the cavity loss k 2 by a factor of (C/A) 2 /2, 
and hence we will neglect it. Similarly the radiative atom damping term is smaller than 
the intrinsic atom decay rate and we will neglect it. In the limit of weak pumping, down 
converted pairs of photons are produced incoherently at a rate 72 2 (C/ K i) 2 given by Fermi's 
golden rule with n 2 the average number of pump photons. 
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To understand the squeezing produced by strong pumping, one needs to consider the 
substantial de-excitation rate 7 from the excited state of solid-state qubits. Furthermore, 
the qubit dephasing rate 7 „ is typicaHy at .east one order of magnitude greater than 7 . For 
charge qubits, optimal phase coherence occurs at the charge degeneracy point [17[ where 
(it happens that) the symmetry prevents the three wave mixing which we require for PDC. 
Hence it is crucial to take into account these effects to obtain a realistic estimate of squeezing. 
For a deviation of 10% from the charge degeneracy point, coherence times of T 2 ~ 200ns 
have been demonstrated We start with the Hamiltonian Hj obtained via the unitary 
transformation UiU 2 in Eq. which is defined in the rotating wave frame of tv 2 , 

a r 

Hj = —a z + 92 sin 9 \ a p \ a y + — (a[a\o~~ + aiaia + ). (8) 

Here T = 2g\ •sm2QjiO\ and we have taken the pump to be classical: (a 2 ) = a p = i\a p \. 

Following the standard quantum theory of damping, we investigate the coupled system 
of qubit and cavity plus the reservoir. After integrating out the reservoir degrees of freedom 
and using the Markov approximation, one obtains the master equation Q for the reduced 
density matrix p of the qubit plus cavity system, 

dp Tr~r 1 r t 1 t 1+171"-+ If- 1 + -1 lv r 1 

— = -i[H I ,p\+K 1 a x pa{ - -a[a x p- -pa[a 1 \+- a pa* - -a ] a p--pa l a +—[o z po z - p\ 

n (9) 

Using the quantum regression theorem[3], the variances V(u) (homodyne spectrum) of 
quadrature X\ = a + a^ and X 2 = (a — aJ)/i for the output cavity photon mode are given by 

V{u) = 1 ± «! J™ c/re-^ T Tr ± a\)e Cr (a lPss ± p ss a{)} , (10) 

where £ is the Liouvillean operator, p ss the density matrix at the steady state, and the 
(+) and (— ) sign correspond to X x and X 2 quadratures, respectively. We have numerically 
calculated the V(u) of quadrature X\ and X 2 based on an exact diagonalization study 
in the Hilbert space of size 2N x 2N corresponding to the two possible spin states and 
the number of photons being restricted to less than N . We have chosen the following 
conservative set of experimental parameters: uji/2it = 3GHz, 7\ = 2ps,T 2 = 100ns, Q = 
5 x 10 5 , g 2 /27r = 18MHz, E c /Ej = 1, and sin 6 = 1 near the CDP, where T 2 = 7J 1 , 
7_i_ = 7^ + 7/2, and A/ui = 0.05[11(. In Fig. |21 we have plotted V(u) as a function of 
uj for A = 2.6 x 10 4 , T = 5, iV = 10, and g 2 \a p \ = 2827 in units of k x = Ui/Q. We find 
that the maximum squeezing is obtained at uj = as expected. The frequency width of 



the squeezing spectrum is controlled by the cavity width k%. The variable V is tunable by 
varying the qubit gate voltage. It is experimentally observed that the dephasing time T 2 
decreases very rapidly, as the qubit is detuned from the CDP Hence we will restrict 
the relative deviations from the CDP to be small |2iV 9 — 1| < 0.1, where T 2 > 100ns. This 
restriction yields a maximum value of T = (8g%/u 1 )(E c /Ej)(l-2N g ) ~ 0.53MHz. In Fig.H 
V(lu = 0) is plotted as a function r for iV = 10, 12, 15, 20. We have checked that AX 1 ■ AX 2 
satisfies the minimum uncertainty bounds and closely approaches the minimum uncertainty 
condition. The maximum output squeezing is obtained by extrapolation of the numerical 
results for finite N to N = 00 yielding V(u = 0) = 0.05, that is, — 13dB, as shown in the 
inset of Fig. El 

We have performed extensive simulations by varying the sets of parameters. Increasing 
A consistently improves the squeezing because it reduces the effect of dephasing of the 
qubit. We can analytically study the effect of finite qubit decay time and dephasing on 
the squeezing as follows. Based on perturbative analysis, we have demonstrated that the 
PDC rate is given by {g 2 T said / A)a z as shown in Eq. When concerned with the spin 
dynamics of the qubit alone, we may neglect the coupling term between cavity photon 
and the qubit, which is much smaller than the other terms in Eq. |HJ By neglecting the 
coupling, we obtain the following steady state solution for the spin polarization (cr z ): (<j z ) = 
~- l(l± + A 2 )/(4g 2 |a p | 2 s i n2 ^7-L + 7(7± + A 2 ))- I n the absence of pump photons, (\a p \ = 0), 
the qubit correctly decays down to the ground state, that is (a z ) = —1. For large detuning, 
the finite pump only slightly quenches the spin polarization. 

Due to the quenching of qubit spin, the effective PDC parameter \2 is given by the 
following formula 



X2 ~ A 4s 2 2 K| 2 sin 2 #( 7± / 7 )+7i + A 2 - 1 1 

Maximum squeezing is achieved near the critical point \2 — K i/2[3j, which leads to the 
optimal value of T for A 3> 7^ 



r* = — 

2x 



l + 4( 7 ±/7)* 2 , (12) 



where x = g 2 sin 8\a p \/A. In Fig. T* is plotted as a function of A. The dotted line 
represents the result from the analytical formula of Eq. ^] which shows excellent agreement 
with the numerical simulation (filled circles). By differentiating with respect to x, one 
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can obtain the minimum value of T* required to achieve the maximum squeezing: r^ in = 
2fti(7_i_/7) 1//2 for x* = (7/47_i_) 1 / 2 . Since the value of the T/ki can reach as large as 13.3 
for the chosen set of parameters, it is well above T^^/ki = 2(7j_/7) 1 / 2 = 8.9. Hence the 
maximum squeezing can be realized with the chosen experimental parameters. The validity 
of perturbation theory in g 2 Q| imposes the following constraint on the product of the 
pump amplitude and the coupling, ^l^pl < (A/2 sin 0). We note that one can control the 
values of g\ and g 2 independently by shifting the position of the qubit within the cavity. 
When the CPB is located at 1/4 of the resonator length from the left edge, g 2 vanishes as 
shown in Fig. ^ Hence by placing the CPB slightly left of the above position, the pump 
amplitude can always be made large enough for the classical approximation to be valid. 

In summary, we have studied degenerate parametric down conversion and squeezing in 
circuit QED, where a superconducting Cooper pair box (CPB) qubit is located inside a 
transmission line resonator. We have shown that away from the charge degeneracy point, 
the coherent coupling of different cavity modes through the qubit can generate a squeezed 
state via three-wave mixing. We have investigated the effect of the finite qubit lifetime 
and dephasing on squeezing, which will be crucial especially for the qubit away from the 
charge degeneracy point. By performing both the numerical and analytical calculations, we 
have demonstrated that the squeezing of about 13dB below the vacuum noise level can be 
obtained for the currently available experimental parameters. 
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FIG. 1: (Color online) The schematic diagram for the circuit QED coupled to a Cooper-pair box 
located at one edge of the cavity. The red and blue lines represent the fundamental and the second 
harmonic cavity modes, respectively. The change of CPB position from the edge will vary the 
coupling strengths g\ and g2- 
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FIG. 2: The quadrature variance V(co) for X\ is plotted as a function of lo for N = 10, A = 
2.6 x 10 4 ,T = 5,g2\a p \ = 2827 in units of k\. The maximum squeezing is obtained at uj = 0. 
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FIG. 3: The quadrature variance V(u> = 0) for X\ is plotted as a function of V for several different 
values of N = 10, 12, 15, 20. The minimum value of V{uj = 0) decreases with the increase of the 
maximum photon number N. The critical value of T is about 8.2 in units of k\. In the inset, the 
maximum output squeezing for several values of N is plotted with respect to 1/N. By extrapolation 
to N = do, we obtain the V m i n (u = 0) = 0.05, that is, — 13dB. 
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FIG. 4: The critical value T* is plotted as a function of A. The dotted line represents the result 
from the analytical formula and the filled circles from the numerical simulation. Here T* and A 
are in units of n±. 
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